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Abstract We construct a dense packing of regular tetrahedra, with packing density D > 
.7786157. 

Keywords Crystallography • Packing • Regular solid • Hilbert problem 



1 Introduction 

"How can one arrange most densely in space an infinite number of equal solids of given 
form, e.g., spheres with given radii or regular tetrahedra with given edges (or in prescribed 
position), that is, how can one so fit them together that the ratio of the filled to the unfilled 
space may be as great as possible?" (excerpt from David Hilbert's 18th problem [8]). 

This is a very old problem. Aristotle [1] believed that you could tile space with reg- 
ular tetrahedra. Everyone believed him for the next 1800 years, until Johannes Miiller 
(aka Regimontanus) proved everyone wrong (histories by Dirk Struik [13] and Marjorie 
Senechal [12]). 

This is a very understandable mistake, because tetrahedra almost tile space locally, but 
not quite. The cluster E 5 (consisting of 5 tetrahedra joined symmetrically about an edge) 
has total solid angle 5 -2cos _1 i ss 12.309594173408 about the edge, and local density 
D rts .979566380077. The cluster V 20 (consisting of 20 tetrahedra joined symmetrically 
about a vertex) has total solid angle 20 • (-7T + 3 cos" 1 j) w 11.025711968651 about the 
vertex, and local density D w .877398280459. 
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For the sphere, the analogous problem was very challenging. A long time ago, Jo- 
hannes Kepler conjectured that the densest packing is the hexagonal close-packing (HCP), 
with density D = k/\/V& w .740480489693. Carl Friedrich Gauss proved that the HCP 
is the densest lattice packing. Only recently, Thomas Hales and Samuel Ferguson [6, 7] 
proved that the HCP is the densest packing in general. 

Helmut Gromer [5] constructed a lattice packing of the single tetrahedron B 1 with 
density D = || w .367346938775. Douglas Hoylman [9] proved that Gromer's packing 
was the densest lattice packing. 

Andrew Hurley [10] constructed the tetrahelix. I calculated its packing density D = 

yf^^j w .531273435694. 

Hermann Minkowski [11] proved that the densest lattice packing of any convex body 
must satisfy certain constraints. Ulrich Betke & Martin Henk [2, 3] developed an efficient 
computer algorithm to compute the densest lattice packing of any convex body, and they 
applied it to the Archimedean solids. I used Betke & Henk's program to calculate the 
packing density of the convex hull of V 20 , D w .716796401602. 

John Conway & Salvatore Torquato [4] used Betke & Henk's program to calculate the 
packing density of V 20 inscribed in a regular icosahedron, D > .7165598. They slightly 
improved the packing density by wiggling the tetrahedra, D > .717455. They raised the 
question whether the packing density of regular tetrahedra is, not only less than that of the 
sphere, but perhaps even the lowest of any convex body. 

In this paper, we will construct a dense packing of regular tetrahedra with packing 
density D w .778615700855, which beats the densest sphere packing! 



2 Description by Construction 

Feel free to look at the pictures in Sec. 3, as you read the descriptions in Sec. 2. 

We construct a 2-parameter family of clusters of regular tetrahedra. Each cluster (9 
tetrahedra) is the union of 1 'central' tetrahedron, 4 'upper' tetrahedra attached to an edge 
of the central tetrahedron, and 4 'lower' tetrahedra attached to the opposite edge of the 
central tetrahedron. The 4 upper tetrahedra can rotate about its attached edge, by an angle 
parametrized by u, and the 4 lower tetrahedra can rotate about its attached edge, by an 
angle parametrized by v. Clusters have 2 orientations: 'positive' and 'negative' (scalar 
mult, by —1). 

For each parameter value (u, v), we construct an optimized packing of clusters, which 
is crystallographic. The clusters pack in layers with alternating orientations. The full 
symmetry group of the packing acts transitively on all clusters. The direct symmetry 
group, which is also the translation symmetry group, acts transitively on a coset. There are 
2 cosets, which correspond to the 2 orientations = 2 translation classes. So the fundamental 
domain (under the lattice of translations) contains 2 clusters =18 tetrahedra. 

Since every cluster is equivalent to all other clusters, we restrict our attention to the 
cluster at the origin and its immediate neighbors. Also, since the cluster (9 tetrahedra) 
is extremely non-convex, we think of it as the union of (the convex hull of) the 'upper 
half-cluster' (5 tetrahedra) and (the convex hull of) the 'lower half -cluster' (5 tetrahedra). 
In order to check that clusters don't overlap in the packing, it's convenient to check that 
half-clusters don't overlap. 



Each cluster has transverse edge-to-edge intersections with neighboring clusters in the 
same layer, and partial face-to-face intersections with neighboring clusters in adjacent 
layers. For intersecting edges, the separating plane contains the union of the edges, and 
for non-intersecting edges, the separating plane is between and parallel to both edges. 
For intersecting faces, the separating plane contains the union of the faces, and for non- 
intersecting faces, the separating plane is between and parallel to both faces. 

The intersections determine equations in terms of the cluster coordinates and lattice 
vectors, which we solve in terms of (u, v) and optimize over all packings in the family. The 
maximum occurs at (u,v) w (-.034789016702, + .089604971413), with packing density 
£>« .778615700855. 

2.1 The cluster B^ v [See Fig. 1] 

We construct the cluster B 1 , consisting of 1 tetrahedron. It is the convex hull of the 'upper 
edge' + and 'lower edge' + 

We attach 4 additional tetrahedra face-to-face along the 'upper edge', and we call this 
cluster E^+. There are now 5 tetrahedra along the 'upper edge', and they form the 'upper 
half-cluster' = B 1 UE^+. 

We attach 4 additional tetrahedra face-to-face along the 'lower edge', and we call this 
cluster E^~. There are now 5 tetrahedra along the 'lower edge', and they form the 'lower 
half-cluster' E* - = B 1 LIE*". 

There are 9 total tetrahedra, and they form the 'cluster' BjJ v = E^ + UE^. Clusters have 
2 orientations: 'positive' clusters +B' and 'negative' clusters — B^ (point reflection). 
The set of clusters {B^ v : — g < « < + 5 , — 5 < v < + 5 } forms a 2-parameter family of 
clusters. 

2.2 The swivel parameters (u,v) [See Fig. 2] 

The cluster E^ + can rotate about the 'upper edge' through small angles, without overlap- 
ping with B 1 . The 'upper swivel parameter' —\<u<+^ corresponds to the x-coordinate 
of the highest vertex g+, and the rotation angle is given by u = \/3 sin 9 U . 

The cluster E^~ can rotate about the 'lower edge' through small angles, without over- 
lapping with B 1 . The 'lower swivel parameter' — g < v < +g corresponds to the x- 
coordinate of the lowest vertex q~, and the rotation angle is given by v = \/3 sin 9 V . 

2.3 The optimized packing P„ v [See Fig. 3] 

The packing is made in layers. The 'central' layer consists of translates of 'positive' 
clusters, in a 2-dimensional lattice \? u v = Z • 2a + Z • 2b with basis {2a, 2b}. (When 
(u,v) — (0,0), this lattice is a square lattice. For our small values of (u,v), this lattice 
is almost a square lattice.) 

The layers alternate between 'positive' layers containing 'positive' clusters, and 'neg- 
ative' layers containing 'negative' clusters. The layer 'above' the central layer consists of 
translates of 'negative' clusters, in the offset c — L% v . The layer 'below' the central layer 
consists of translates of 'negative' clusters, in the offset d — L% v . 



The 'positive' clusters live in a 3 -dimensional lattice v = v + Z(c — d) with basis 
{2a,2/?,c — d}, and the 'negative' clusters live in the offset c — L^ v = d — L^ v . The 
packing consists of 2 cosets of the lattice v , which correspond to the 2 orientations. For 
each cluster v , we find the optimized packing and call it P !( v . 

The packing is periodic and crystallographic, and has a large isometry group which acts 
transitively on clusters. There is a direct isometry between any two clusters of the same 
orientation, and an indirect isometry between any two clusters of opposite orientations. 
The direct isometry group is the translation group Lf ( v = Z ■ 2a + Z • 2b + Z(c — d). The 
point group is (generically) trivial. 

The set of optimized packings {P u v : — ^ < m < + g , — \ < v < + g } forms a 2-parameter 
family of packings. 

2.4 The neighbors & separating planes of v [See Fig. 4] 

Since B^ v is non-convex, we consider the half-clusters Ef ( + ,E^~ separately. Since each 
half-cluster is disjoint from its neighbors, the union is also disjoint from its neighbors. 
Therefore, we can verify that we have a packing. 

Separating planes separate the upper half-cluster Ef, + from each of its neighbors: the 
plane R • 2a + R • 2b + for the upper above half-layer, 4 planes for the 6 neighbors in the 
lower above half-layer, 4 planes for the 8 neighbors in the upper central half-layer, 6 planes 
for the 8 neighbors in the lower central half-layer, and the plane R • 2a + R ■ 2b + d + 
for the upper below half-layer. (g+ is the highest point of the upper half-cluster.) 

Separating planes separate the lower half-cluster from each of its neighbors: the 
plane R ■ 2a + R • 2b + c + q~ for the lower above half-layer, 6 planes for the 8 neighbors in 
the upper central half-layer, 4 planes for the 8 neighbors in the lower central half-layer, 4 
planes for the 6 neighbors in the upper below half-layer, and the plane R ■ 2a + R • 2b + q~ 
for the lower below half-layer. (q~ is the lowest point of the lower half-cluster.) 

2.5 The intersections & separating planes of B„ v [See Fig. 5] 

In the same layer, clusters touch transversely edge-to-edge, intersecting at only 1 point. 
Each upper half-cluster has 4 such contacts with neighboring upper half-clusters, and each 
lower half-cluster has 4 such contacts with neighboring lower half-clusters. These inter- 
sections give 4 equations. Each upper half-cluster has < 2 such contacts with neighboring 
lower half-clusters, and each lower half-cluster has < 2 such contacts with neighboring 
upper half-clusters. These intersections give < 2 equations. 

For non-intersecting skew edges E[a,b] ^ E[c,d] + w, there are many choices for the 
separating plane. We can choose the plane which contains E[a,b] and parallel to E[c,d) 
(ie: the plane through the points a,b,b + c — d), the plane which contains E[c,d] and 
parallel to E[a,b] (ie: the plane through the points c,d,d + a — b), or any parallel plane in 
between. 

For the transverse edge-to-edge intersection E [a , b] = E [c, d] + w, the separating plane 
contains the union of the edges E[a,b] UE[c,d]. We solve the virtual intersection equation 
sa + ( 1 — s)b = tc+ (1 — t)d + w, between lines with parameters s, t . For our (small) values 
of (u, v), our edges intersect in their interiors (ie: < s < 1 and < t < 1). 



In adjacent layers, clusters touch partially face-to-face. Each upper half-cluster has 3 
such contacts with neighboring half-clusters in the above layer, and each lower half-cluster 
has 3 such contacts with neighboring half-clusters in the below layer. These intersections 
give 6 equations. 

For non-intersecting parallel faces F[a,b,c] ^F[d,e,f] +w, there are many choices for 
the separating plane. We can choose the plane containing F[a,b,c], the plane containing 
F[d,e,f], or any parallel plane in between. 

For the partial face-to-face intersection F[a,b,c] — F [d,e,f] + w, the separating plane 
is the plane which contains the union of the faces F [a,b,c] U F[d,e,f\, There are two 
ways to write the intersection equation. 

We can reduce it to a transverse edge-to-edge intersection of any non-parallel pair of 
edges from {E[b, c] , E [c, a] , E [a, b] } and {E [e, f] , E [/, d] , E [d, e] }). If we choose carefully, 
we can find two edges which intersect in their interiors. 

Alternatively, we can reduce it to an incidence between a triangle F[a,b,c] and a vertex 
from {d,e,f}, or vice versa. Then we solve the virtual intersection equation sa + tb + 
(l—s — t)c = d + w, with plane parameters s, t . If we choose carefully, we can find a point 
which lies inside the triangle (ie: < s < 1 and < t < 1). 

3 Pictures & Equations 

The figures are really super-figures which contain multiple sub-figures. Figs. 1-5 are 
floor-plan diagrams which show how the various layers fit together (ie: the top layer 
is placed at the top, and the bottom layer is placed at the bottom). For the purpose of 
perspective, we draw the coordinate axes as solid objects. For the purpose of comparison, 
we provide both top views (perspective vector (0,0, 1)) & side views (perspective vector 
(0, -cos ^, sin ^)). 

• Fig. shows the orientation of the (x,y,z) coordinate axes (for Figs. 1-5). 

• Fig. 1 (2 pages) shows the cluster Bl v . 

• Fig. 2 (2 pages) shows the swivel parameters (m, v). 

• Fig. 3 (2 pages) shows the packing P M v . 

• Fig. 4 (2 pages) shows the neighbors & separating planes (general overview). 

• Fig. 5 (4 pages) shows the intersections & separating planes (fine details). 



Fig. (I 1 The (i.Vi) coordinate 
axes {side view & top vi&v). 
for Figs. 1-5 




B 1 : 1 tetrahedron 

coordinates: fi[<+ 1 ,+ 1 .+ 1 ),(+ 1 - 1 ,- 1 >.(- 1 ,+ 1 - ] >,<- l.-l.+l >] 



Clusters have 2 possible orientations: 

positive' or 'negative' (point reflection = scalar mult, by — 1 ) 



E„ : 4 tetrahedra, joined face<-»face 

attached to the 'upper edge' £[<+! ,+ I 1 ,- 1 ,+ l)| 

coordinates: S[<+ 1 ,+ 1 ,+ 1 ),(- 1,-1,+] )p*J>*] 
8t<+l,+ l,+ l>.(-l.-l.+ l>^v?;] 

s[<+ 1, + 1. +i ),<-!,- 1, +i>^„v:] 

B[<+1,+ 1, + !),{- 1,-1,+ 

K = - 2 „~ + l„+i ^ T77, + 1 + i u- 7 - 4~6 - 2 » 2 ) 

pi = <+ V 6 - 2 m 2 6 - 2 ir ,+ 1+f «+} V 6 - 2 m 2 ) 

f/+ = (+H.-M.+ 1+^6- 2 u 2 ) 

,■+ = (+^+f V 6 - 2h 2 ,-ju-jV 6 " 2 " 2 ,+ 1-fw+yV 6 " 2 " 2 > 
4 = <~«+|V 6 - 2 « 2 ,+iu- 4 -^6- 2u* ,+l-|„-lV6- 2 » 2 > 




Fig. 1^ The positive cluster +B^ r w"«H>) 



E^~ : 4 tetrahedra, joined face<-»face 

attached to the 'lower edge' £[<+! -I 1 ,+ 1 ,-l>] 

coordinates: fi[<+ 1,-1 ,-1)^-1 ,+1,-1 )jo~,p~] 
fi[<+ 1, - 1,-DX- 1, + 1,-1>,/;-^] 
B[<+ 1,-1, -!),{- 1, + 1,-1)^^] 
B[(+ 1 ,- 1 ,- 1 ),{- 1 ,+ 1 ,- 1 )r~ j~] 

o- = <-Iv,-iy 6 - 2 v> , ->-|V & - 2 ^ . -1-^+^ 6 - 2 v 2 ) 
/j; = <+ fv-f V 6 - 2 v 1 ,+}v-f V 6 " 2 v 2 v 4 V 6 " 2 ^ > 

q~ = (+v.+v,-l-V6-2v 2 ) 

,■- = (+ ji'+f V 6 - 2 v 2 6 - 2 ,-1+yV-j-y 6 - 2 > 

s; = <-lv+i V 6 - 2 v 2 ,-ly+f V 6 - 2 ^ ,- 1+f v+| V 6 - 2 v 2 I 




Fig. 1^ The negiitivt cluster - B* y (yjrfe v&w) 



/i 



Inside the cluster E^ + , there is a small gap 

around the 'upper edge' E[(+ 1 ,+ l ,+ 1 ),(- 1 - 1 

-> the cluster can rotate about B 1 , by small angles 

11 parametrizes the (x,— y) ■coordinates of q* (indirectly, the angle) 
-> the endpoints of 11 are determined by intersections 

(+1 -i -i>-» 




Fig. 2siik- The swivel parameters (u.v) (side view) 



Inside the cluster E^ - , there is a small gap 

around the 'lower edge' E[(+ 1 ,— 1 ,— 1 ),(— I ,+ 1,— 1)J 

-» the cluster E^!~ can rotate about B 1 , by small angles 

V parametrizes the (x,+y)<»ord i n ates of (indirectly, the angle) 
-» the endpoints of v are determined by intersections 

j- = (+l,+ l,+ l>-, v = +| 




Fig* 2tcvp The swivel puriimeters {u,v) (top view) 



c 



Clusters pack very efficiently in layers 

the ZD-lattice L^ v has almost-square basis {2a,lb} 



Each layer is identical to any other layer (under some isometry) 
-> we need to examine only the layer at the origin, and its neighbors 

Z > : 'above layer' at offset c 

Z - : 'central layer' at offset 

Z < : 'below layer' at offset d 



above 



central 



below upper below 




Fig. 3,i,t. The packing P„. r {side view) 



-> the 3D-lattice has basis {2a,2b J c—d\ 



Layers have 2 possible orientations: 

'positive' layers contain "positive' clusters 
'negative' layers contain "negative' clusters 




Planes separate 
the half-cluster E* 4 from 
all of its neighboring 
half-clusters E' + & 
-> there is no overlap 



-> planes separate 
the cluster B* from 
all of its neighboring 
clusters B' t . 
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plane : (2a^b}+q* 



upper above 




6 neighbors 




Fig. 4* The neighbors & sepiirating planes of the upper half eluster E;j + 




1 

lower above 




The faces between 
the 'upper central' and 
the 'lower above' layers 
have opposite orientation 
(point reflection = 
scalar mult, by -1) 
-» the faces are parallel 
-» the intersections are 
partial face <-> face 

v -2a 
*-2h 



sep «*•>. 

sep E*V- 
sep El + r- 

5+ 



piane thru F\ (+ 1 ,+ 1 .+ 1 ) #* \#-f\{+\ ,+ 1 ,+ 1>^I /*] 
plane thru F\ <- 1 .- 1 .+ 1 > //* j% | x-F\ ( - 1 ,- 1 .+ 1 > 4* ,p+ 1 -2a 
plane thru F\ <- 1 .- 1 .+ 1 > & f* \,c-F\ ( - 1 ,- 1 ,+ 1 >^/J |-2A 

planethruF|{+[,+ l.+ IXX| II c-F\(+ I > +l,+l)j^*,p*]-2n+2ft Sep E* + ^-Ej + -2a+2ii 
plane containing (the union of) intersecting faces is well-defined 
plane parallel to parallel faces has 1 degree of freedom between the faces 

Each cluster is congruent 
to any other cluster 
(under some isometry) 
-* we need to examine 
only the hall cluster E;j + 
22 neighbors 
14 separating planes 
9 intersection points 
7 equations 

K'la-I, plane thru E[{+\ ,+ 1 ,+ ]>/+].£[<- 1 ,- I ,+ l).p+]+2a-2t> sep E;|\E?, + +2«-2i> 
K"_ 1Hl , plane thru £[<-] .- 1 .+ ] ),p^\£[{+ 1 .+ 1 .+ 1) /■* \-2a+2b sep El+£l + -2a+2b 
Kjj plane thru E[(- 1 .- 1 .+ !)/£],£[{+ 1 ,+ 1 ,+ 1 >./>,; |-2J» sep E^ .Ef, + -26 
K° A plane thru £«+],+ ] ,+ l>^J]JE[<-l,-l,+l>/+]+ia sep E* + ,Ef, + +2fc 
plane containing (the union of) intersecting edges is well-defined 
plane parallel to skew edges has 1 degree of freedom between the edges 





plane thru E[(+],+ ].+ Y}j>*]£[p;,o;.\+2a 
plane thru E[(-\ .-] .+ ])#+) || £[c,:jv]-2o 
plane thru E\r* ,s*\ || E[(- 1 .+ 1 .- I>^~]— 2A 
plane thru ,£>*],£[(+ 1 -1 -\),o;,]+2b 



L'Vj, Plane thru /=■[<- 1. - 1. + 1 ).<- 1, + 1 ,-l),< + 1.-1.-1)] 
L+„+j, plane thru F[<+l.+ l.+ l),<+I.-l -1 ).<-], + 1,-1)] 
Mr. The i merseciions & separating planes of E** {side view) 



sep Er^f 
sep 

sep E^.Ej- 
sep E^" 



+2a 
-2a 
-2b 
+2b 



sep E;|\Ef -2«-2i> 



sep E^^' 



+2a+26 




GJ fl <+ 1 ,+ 1 ,+ 1 > 4l< rl j = c-F[(+ 1 ,+ 1 ,+ 1 > 4+ j-U 

g* f[(-i,-\,+\)x/:. p :\ = c-F[(-\.-\.+\)4„ P :\-2a 

c; f[ (- 1 .- 1 .+ iVaj/j 1 ■ c-F[(-i .- 1 .+ 1 >^ /J]-2i> 

G, + A fl (+ 1 .+ 1 .+ 1 ) 4: ] * c-F[( +l.+ l.+ Dvg:^]-2a+2* 
3 intersections (partial face «-> face) 
1 non-intersection (parallel faces) 





G»_„ £[(+ 1 ,+ 1 ,+ 1}j-*| = £[{-1 .-] ,+ \)tf]+2a-2b 

gju £[<- 1 .- 1 -+ >>,p,;i = £[<+ !.+"■+ i)K\-^+2b 

G'l E[<- 1 .- 1 ,+ !>/•♦] = £[<+ 1 .+ 1 .+ 1 ),pl\-2b 
G° £[(+l,+J,+ L)^] = E[<-l,-I,+ I)/-;]+2* 
4 intersections (transverse edge*->edge) 
non- intersect ions (skew edges) 



2 interseetfi 




Hi £[<+l, + !,+ !),<] - E[p;»;]+2a 
H; Et<-l~l,+l>^] * E[r-j;]-2a 
Ht E[r* ] * E[(- 1 .+ 1 .- 1 ) j~]-2b 

n: h e[ P : ,0:] =£i<+i,-i,-i> A :]+2& 

2 intersections (transverse edge<->edge) 
2 non-intersections (skew edges) 

Fig. 5* The intersections & separating planes of E^ r {top view) 



o 





/ 6 separating planes 




I + 

K 

v: 



I, 

a i b 



plane thnj £[<-l ,+ 1 ,- 1 > T s;] || £[i-+^]+2fi 
plane thru £[<+] ,-] -]),o;],E[p*,o*,\-2b 
plane thai g[r;^] || £l(-l ,- 1 .+ \)jt>Z\+2a 
plane thru E[p~ /)~\£[(+\,+ \ ,+ l)^j-2o 
plane thru £[<+! -1 -1 >,<+! .+ 1 ,+ 1>,<- 1 -1 .+ 1 >] 
plane thru £[{- 1 ,+ 1 ,- 1 ),{- 1 ,- 1 ,+ 1 },{ + 1 ,+ 1 .+ 1 )J 




sep E^,Ej + +2f> 
sep Ef",Ef, + -2/; 
sep Ef",Ef,*+2« 
sep E^j:f, + -2fl 
sep Ef-,E; + +2w-2ft 
sep E*"Sj + -2a+2ip 
Each cluster is congruent 
to any other cluster 
(under some isometry) 
-* we need to examine 
only the hall-cluster E;*~ 
22 neighbors 
14 separating planes 
9 intersection points 
7 equations 
sep E*" 



,E" i-2a + 2h 
sepE^JE?--2«-2/j 
sep E'" J&l~+2a 
sep E*",lj"-2a 



K^t, plane thru £[<-l,+ l,-l)^],£K+l ,-1 ,-l), P ;]+2«+2fc 
K-a-fc plane thai £[(+1 .- 1 .-1 ), P ;\£[{- 1 .+ 1 .- 1) j~; \-2a-2b 
K° +tl plane thru £[<+ 1 ,-1 ,-\)s;),E\<- I .+ I ,- I >./;" |+2« 
plane thru £(<-] ,+ 1 ,-\),p-]£[(+l ,-1 r;\-2a 
plane containing (the union of) intersecting edges is well-defined 
plane parallel to skew edges has 1 degree of freedom between the edges 

The faces between 
the 'lower central' and 
the 'upper below' layers 
have opposite orientation 
(point reflection - 
scalar mult, by -1) 
-» the faces are parallel 
-» the intersections are 
partial face*-»face 
sep E'^^-E*" 




K5 
K 



plane thru F\ (- ] .t 1 ,- 1 ) ff, ?' Ut~F\ ( - 1 ,+ ] ,- 1 ) jy" j;| 
plane thru F\ {+ 1 .- 1 .- 1 ) x/; ,p~ \ M-FH + 1 .- 1 .- 1 .;>; \~2h 
plane tbni .- 1 .- 1> s~]ji-F[(+l - 1 .- 1 )jq; f~\+2a 
plane thru fl<-[ .+ 1- l>^T.p;| || rf-f [<-!.+ 1 -I ]-2a-2b Sep E^," jcl-tf.~-2u-2b 



sepE--.f/-Er-2t 
sep Ef-^-E^+2fl 



plane containing (the union of) intersecting faces is well-defined 

plane parallel to parallel faces has 1 degree of freedom between the faces 

Kg. ?"„;.. The i iiUTseetions & separating planes of E'~ (side view) 




h: 



£[<-!.+ l.-l>j-J *E[r + u 4\+2h 
Et<+l,-J,-l>^]«£|»^-2i 
E[r;,s;\ * E[(-[-lM),o+]+2a 
E[p^ I = E[{+ ! ,+ 1 ,+ 1 )j*]~'2m 
2 intersections (transverse edge<->edge) 
2 non -intersect ions (skew edges} 



oh 





























4 intersections 






G° a+h £[<-!,+ 1 ,-1>j--J = £[(+1 .-1 .-\),p;\+2 a +2b 
G° a+h £[(+ 1 .- 1 .- 1 }. P ;] = £[<- 1 ,+ 1 .- I ).r;.\-2a-2b 
Gj £| <+ 1 .- 1 . - 1 } f - 1 = £[( - 1 .+ 1 .- 1 ),p- l+2a 
G« £l<- 1 ,+ 1 1)^-] = £l<+ 1,-1-1) ,r;]-2a 
4 intersections (transverse edge*-»edge) 
non -intersect ions (skew edges) 





£[(-l,+ l -1)^/-;]: 

FK+i.-i.-i)^: ./>,-] 



-.d-F\(+i.-i-[).q;, P ;.)-2b 



««, 

G; £| (+ 1 .- 1 .- 1>^~ r~ | = d-F[ (+ 1 .- 1 .- 1 )4~ /fl+2a 
G; h F[(-i,+ L-[^;,p;]*d-F[{-i,+i-i)^;.,p;}-2a-2b 

3 intersections (partial face*-^face) 
1 non-intersection (parallel faces) 

Fig. 5,V, The intersections & separating planes of Uop view) 



n 



4 Proof by Computation 



4. 1 Numerical results 

We use the following procedure to compute the maximal optimized packing: 

1. Write abstract expressions for the virtual cluster vertices {o+ ,q„ ,r+ ,s„}, 
{o,7 , p~ , q ~ , r~ , s~ }, in terms of the swivel parameters (u,v). [see Fig. 1] 

2. Write abstract equations for the virtual intersections between the cluster and its neigh- 
bors, in terms of the virtual cluster vertices (ultimately in terms of the swivel parameters), 
the lattice vectors {a,b,c,d} and extra intersection parameters (edge-to-edge or point-in- 
plane), [see Fig. 5] 

G = { G ' G ai G t> G ' G a ' G fo > G a> G b' G a+b' G a-b}- 
H = {H a ,H fc } 

3. Solve for the abstract lattice vectors and extra intersection parameters, in terms of the 
swivel parameters. Clear the common denominator. 

4. Write abstract expressions for the lattice volume V & packing density D, in terms of the 
lattice vectors (ultimately in terms of the swivel parameters). (The fundamental domain 
contains 2 clusters =18 tetrahedra. Each tetrahedron has edge length \/8 and volume |.) 

V = det[2a,2b,c-d] 
D = 48/V 

5. Find the relative minimum of the lattice volume, with respect to the swivel parameters. 
(We could find only a numerical approximation, but not the minimal polynomials, because 
the degree is too big for our computer program.) 

(u,v) w (-0.034789016702, +0.089604971413) 

6 + . Evaluate the vertices of the upper half-cluster. 

o+ w (-1.061384364770, + 1.061384364770, -0.935697925928) 
p+ « (-1.644260072209, + 1.644260072209, + 1.769946511051) 
q+ « (-0.034789016702, +0.034789016702, +3.448995599962) 
r+ w (+1.621067394407, -1.621067394407, + 1.862717222257) 
s+w (+1.115500612974,-1.115500612974,-0.873850785124) 

6~. Evaluate the vertices of the lower half-cluster. 

o~ w (-1.156897066820, -1.156897066820, +0.822958681256) 
p- w (-1.600938143111,-1.600938143111,-1.934876528675) 
q~ w (+0.089604971413, +0.089604971413, -3.446209700372) 
r- w (+1.660674790719, +1.660674790719, -1.695929938240) 
s- w (+1.017511555733, +1.017511555733, +0.982256408212) 



7. Evaluate the extra intersection parameters and intersection points, in order to verify 
that the intersections are valid (i.e.: parameter values between and 1). 

8. Evaluate the lattice vectors. 

a w (+1.286101228632, +0.200477756509, +0.1 17304750804) 
b w (-0.216654244567, + 1.293299854677, +0.064005626691) 
cw (+1.049828831839, +0.572626599359, +4.643104387948) 
d w (-0.331247602500, +1.312711170924, -4.452573714586) 

9. Evaluate the lattice volume & packing density. 

Vw 61.647870634123 
D w 0.778615700855 

4.2 The actual density function D [see Fig. 7] 

For each cluster v , we compute the optimized packing P u v and packing density D M v . 
The set of all packing densities {D u :— 5<m<+^, — ^<v<+g} forms the 'actual' 
density function D. The range of the graphs is +.7767 < w < +.7787 for Fig. 7 and 
+.7675 < w < +.78 for Fig. 8. 

All packings in the family satisfy the set of equations G, which gives all packings the 
same geometric structure and lattice type. (For small values of the swivel parameters 
(u, v), the intersections from the equations G remain valid.) 

However, when we try to find the maximal optimized packing using only the equations 
G, there are 2 places where the clusters overlap. In this case, we need to impose the 2 
additional equations H in order to prevent overlap. 

For a general cluster v , only a subset of the equations H is necessary. So there are 4 
possible optimized packings using the equations G, GU {H a }, GU{HJ, GU{H fl ,Hj, 
which gives us 4 possible packing densities v , DJJ v , D* v , D^ fo v . 

A priori, we don't know which set of equations to use for a given cluster, or which 
clusters satisfy a given set of equations. So we naively assume that all 4 possibilities are 
valid for all clusters, which gives us 4 'virtual' density functions D°, D fl , T> h , D ab . 

4.3 The virtual density functions D°,D a ,D f, ,D afo [see Fig. 8] 

For each swivel parameter value (u,v), we find the 'valid' virtual density function (the 
set of equations which gives the densest non-overlapping packing for the given cluster). 
Conversely, for each virtual density function, we find the 'valid' region of swivel param- 
eter values (the set of clusters having densest non-overlapping packings which satsify the 
given equations). 

The 4 virtual density functions D° ,D" ,D b ,T> ab are smooth, and their surfaces are pair- 
wise tangent along their intersection curves D° n D a , D° n , D fl n T) ah , T> b n T> ab . All 
4 surfaces are mutually tangent at their common intersection point D° n D a n T) b n D ab ', 
which is also the common intersection point of the 4 intersection curves. 



Fig. ft The: (b,v,iv) coordinate 
axes {side mw & lop view), 
far Figs. 7-8 




Fig. 7 The composite density 
function D with {u,v) -traces 
& u''Contours 




The actual density function D consists of the union of the 4 valid regions of the virtual 
density functions, which are bounded by the intersection curves. The 4 valid regions 
are pairwise tangent at the intersection curves, and mutually tangent at their common 
intersection point. Thus, the actual density function is also smooth. 

The common intersection point occurs at 

(u,v,w) « (-0.037320921073, +0.033926596665, +0.778365087767). 
The relative maximum of the actual density function D occurs at 



(K,V,W> 



(-0.034789016702, +0.089604971413, +0.778615700855) 



• Figure 6 shows the orientation of the (u,v,w) coordinate axes (for Figs. 7-8). 

• Figure 7 shows the actual density function D. 

• Figure 8 (2 pages) shows the virtual density functions D° ,D" ,D b ,D ab . 



4.4 The highly symmetric cluster and packing P sym 
The highly symmetric cluster Bq has symmetry group l xyz . 



{I,P,P 2 ,P 3 .Q,QP,QP 2 ,QP 3 },P 



" -1 ~ 




' + 10 " 


+ 10 


,Q = 


0-10 
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The highly symmetric packing P sym also has symmetry group I vyr The lattice has a 
square basis, with lattice vectors of the form 



a = 
b = 
c = 
d = 



(+2/,+2./,0) 
(-2/,+2i,0) 
(+i,+j,+k) 
(-j,+i,-k), 



i = ^(-168 + 106^6) w 1.290787503310 
j= i(-88+42\/6) w 0.209557312632 
k= ^(-262 + 238^) w 4.520824771583 



and lattice volume & packing density 



y= 1 

v 357911 

D = 



(-730200320 + 307139840V6) w 61.846569901642 
1 ( 17 11407 + 7 19859 V6) « .776114181859. 



4477040 

The highly symmetric packing P sym is nof a member of the family {P u ,,}, because adja- 
cent layers have only 2 face-to-face intersections. 

The optimized packing P has trivial symmetry group. As with all members of the 
family, adjacent layers have 3 face-to-face intersections. In order to get this additional 
intersection, we need to shift adjacent layers by a small distance (relative to the lattice), 
which destroys symmetry. 



4.5 Isometries of the cluster B„ v and packing P M v [see Fig. 9] 
We can apply the isometries l X yz m several ways: 

1 . Apply the isometries to the optimized packing P H v (which also applies the isometries 
to the clusters and intersection equations). The 8 clusters I x vz • B^ v are congruent to each 
other, the 8 optimized packings l x ,y,z ' P« v are congruent to each other, and the 8 optimized 
packing densities are equal to each other. 

Each optimized packing corresponds to a point of an actual density function. The 8 
different (but analogous) points have the same w-coordinate (packing density), but each 
point lives in a different (but congruent) surface (actual density functions I M , v ,w ■ D). (The 
isometries I M , v ,w are similar to the isometries l x ,y, z , except that the w-coordinate (packing 
density) is always positive.) 
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Fig <Ci The virtual density functions tP Df J? &* with (w.v)-Caces (side view) 
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\,v,w = {I^tR 2 ,^ ,S,SR,SR 2 ,SR 3 }, R 

2. Apply the isometries to the cluster ,„ but fix the intersection equations G U H. The 8 
clusters I™ z ■ v are identical to 

fT}9 R 9 ij9 R 9 ij9 n 9 r9 r 9 \ 

\ D +u,+vi D -\\ : +ui D -u,-\n D +v-ui D +v,+ui D +u-\n D -v-ui D -u,+vS 

and congruent to each other. The 8 optimized packings 

iPPPPPPPPT 

I* +«,+V' i — V,+«> — U, — V ' +v, — K> +V,+M» +«, — V' 1 — V, — U ' — U,+VJ 

are nof congruent to each other, and the 8 optimized packing densities 

{D+ U:+V ,D_,,,+ U ,D_ M _ v ,D +v _ M ,D +V)+14 ,D +JJi _ v ,D_ V) _ M) D_ B)+v } 

are not equal to each other. 

Each optimized packing corresponds to a point of an actual density function. All 8 
points live in the same surface (actual density function D), but each point has a different 
w-coordinate (packing density). 





Fig 8)1, The virtual density functions D jy'JJ'J)" with iv-contours [top Wew) 



3. Apply the isometries to the intersection equations G U H, but fix the cluster B^ v . This 
is equivalent to fixing the intersection equations and applying the inverse isometries to the 
cluster. Since I r y , is a group, the set of packings & densities is identical. 

You could also apply the isometries individually to each vector quantity in the equa- 
tions. For the virtual cluster vertices {o+ ,q+ ,r+ ,s„}, {o~ ,p~ ,q~ ,r~ ,s~}, use the 
highly symmetric cluster B 00 as a paradigm. For the abstract lattice vectors {a,b,c,d}, 
use the highly symmetric packing P sym as a paradigm. 

In any case, we don't get any additional packings by considering additional isometries 
of the cluster, intersection equations or packings. Thus, it suffices to optimize over a single 
generic cluster B„ v , a single set of equations G U H, and a single actual density function 
D. 

The 8 congruences of the cluster B^ v are symmetric (about (u, v) — (0,0)) as a whole, 
so the 8 sets of intersection equations G U H are symmetric as a whole, and the 8 con- 
gruences of the density function D are symmetric as a whole. However, any single set of 
intersection equations is individually asymmetric, so any single density function is indi- 
vidually asymmetric. 



• Figure 9 (2 pages) shows the isometries of the cluster B„ v and packing P u ,,. 




Opposite face/edge pairs: facc<-»face non-equation G* =s edge<-»edgc equation 
(non-intersecting) "far' fates at top of Ef, + as (intersecting) near' edges at bottom of E|j* 

Go F[{+l,-fiMUl/i]*e-FU*iMM)^^] H,; E\r* .<) = £(<-!. + 1.-1) />," }+2a 

g; Fu-\.-L+v,4;.i>:]* l -F\<-\„-\Mu,:.i>;,]-2a h; £[/>>;] =£i<+i ,-i.-i>^~j— 2« 
g; 6 F[<4i.-M,+ixxi*'-^[<+i-+i.+'i4;.rt]-2 iI +2ft h,;, £[p>;i = £[<+i.-i,-i),o-]+2i> 

Oppt)site face/edge pairs: face<->face non-equation GJT a edge<->edge equation H~ 
(non-intersecting) "far' faces at bottom of Ej" a (intersecting) 'near' edges at top of E;J~ 



G„ F[<-l,+ l,-l>^-]*rf-F[<-],+ 1,-1)^-^1 



H,7 E[ ( ; Ji :|,£K-l,-l,+ l),0+2J 



G;7 Fi{+\.-\.-\Uf;.p;.]*d-F\<+i-\.-\)ij;,p;\-2b hj, £1^;^;] = £[<+i,+ i,+ i),o+]-2t> 

G; ft<+l,-l,-l>^:jvi* i /-F[{+l.-l,-l.).^. r -|+2« H,; £[r"^| = £[<-! ,-l.+ l>.f^J+2a 
G^ F[<-i,+i,-i>4^-l*<f-fl<-L.+ i.-i) 4 ;j»;|-2fl-26 H~ A £[p>;] =£[<+!,+ l,+ l>j+]-2a 




Fig. S* 5 * The 8 sets of intersection equations G IJ H under the isomelries 1, 



4.6 Remarks: Dense Packing 



References 



This construction gives a dense packing, for several reasons: 

1 . The half-cluster E 5 (5 tetrahedra joined face-to-face along a common edge) is very 
dense locally. The local density (fraction of solid angle around the edge) is D w .979566380077. 
The relative density of the half-cluster inside its convex hull is D = ^ w .973557692308. 

2. Opposite edges of a regular tetrahedron are skew perpendicular. Inside the cluster B 9 „, 
the two half-clusters are oriented perpendicularly. So the clusters fit together very nicely 
in layers with an almost-square basis. 

3. Between adjacent layers, the half-clusters and their neighbors are point reflections 
(scalar mult, by —1). Therefore, the faces are parallel, and the intersections are partial 
face-to-face. 

4. For small values of (u,v) s=s (0,0), the clusters have an extra almost-symmetry. Be- 
cause the layers have almost-square basis, point reflection (scalar mult, by —1) is almost- 
equivalent to rotation by | radians about the z-axis. Thus, the layers fit together very well 
to form a dense packing. 

5 Conclusion 
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For the B 9 cluster, if we rotate the 8 non-central tetrahedra independently about the central 
tetrahedron B 1 , we can generalize to an 8 -parameter family B a h c d e ^ h = B 1 U cd U 

Intuitively, the maximal packing occurs when there are no gaps among the 4 upper 
tetrahedra, and no gaps among the 4 lower tetrahedra (ie: the special case of our 2- 
parameter family B 9 U V = H^ uuuvvvv , where the non-central tetrahedra form E^+ = E£+ u u 

and E 4 V ~ = E4- v v ). 

From a geometric perspective, it's easy to see why. If you start with a gapless cluster, 
and try to open a gap somewhere, you will only decrease the packing density. The distance 
between neighbors in the same layer will increase, and the perpendicular distance between 
layers will not decrease. Therefore, generalizing will not help us. 

I found packings for other clusters of tetrahedra B 17 , V 20 , E 15 , B 13 , B 9 with packing 
density D > j^, but the winner is B 9 . 

However, it's possible that other arrangements of tetrahedra yet undiscovered could 
give a denser packing... 
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